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ABSTRACT 

A first-principles method, based on density functional perturbation theory, is presented 
for computing the leading order tunability of high-dielectric-constant materials. 

INTRODUCTION 

High-dielectric-constant materials are of great interest for high-performance electronic 
devices. Besides the high dielectric constants themselves, which could in principle allow 
miniaturization of devices beyond that possible with silicon dioxide, tunability of the di- 
electric response by a dc electrical bias opens up a wide range of additional applications. 
Understanding the fundamental physics of high-dielectric-constant insulators does, how- 
ever, pose some significant challenges. In particular, the effects of homogeneous electric 
fields on the electronic and crystal structure of insulators, and on their dielectric and piezo- 
electric properties, involves subtleties that have only recently begun to be systematically 
addressed. Key advances include the correct formulation of the electric field perturbation 
within density-functional perturbation theory (DFPT) E], and the expression of polariza- 
tion as a well-defined bulk property through the Berry-phase formalism Q. Subsequently, 
the effects of finite fields have been formulated using Taylor expansions around zero field H. 
Using a closely related approach, a formalism in which polarization is the independent vari- 
able has been developed 0] which includes a method for computing the nonlinear dielectric 
susceptibility. 

In this paper, we present a first-principles method for the efficient direct computation 
of the tunability of high-dielectric-constant materials to leading order in the electric field. 
This method is designed for implementation into existing density functional pertubation 
theory packages, such as ABINIT g] and PWSCF/PHONON 0, which typically calculate 
second derivatives of the energy and can be extended to compute third derivatives without 
significant additional computational effort. If third derivatives are available, the collinear 
susceptibility of a ferroelectric material can be computed with the sole approximation that 
the electric-field derivative of the electronic susceptibility be fixed at the value at zero field 
with corresponding equilibrium structural parameters. For noncollinear tunability or tun- 
ability of paraelectrics, or if the DFPT packages are to be used for collinear susceptibility of 
ferroelectrics without extension to computation of third derivatives, nearly all of the neces- 
sary higher derivatives can be obtained from finite differences with a calculation at a single 
additional structure. However, to complete the calculation in these cases, it is necessary to 
estimate the highest derivatives of the electronic susceptibility with respect to the structural 
parameters, as their complete specification would require additional finite-difference calcu- 



lations. These approximations should be reasonable for the materials of interest, in which 
the dielectric constant is dominated by the large lattice contribution. 



FORMALISM 

We consider the dielectric response of a polar insulator to a homogeneous static electric 
field £. The first steps of the analysis are identical to those presented in Ref. ||. Given a 
particular unit cell, let R be the atomic coordinates and rj the homogeneous strain relative 
to a suitable reference structure. The total energy E(£) and equilibrium structure R(£) 
and r)(S) of the system in the field £ are obtained by the minimization 

min£(R,r/,£) (1) 

with R(£), rj(£), and E{£) being the values at the minimum. The polarization P(R, rj,£), 
the thermodynamic conjugate of £, can then be obtained from the expression P(R, rj,£) = 
~ ag 7 '^ an d obtained by evaluating P(R(£), r)(£), £), or equivalently, - dE ^ . 

In DFT, there is as yet no rigorous formulation for finite £, so we expand in £ around £ 
= 0: 



+ ^£ a £pE^(R, rj) + -£ a £ £,E^(R, v ) + ■ ■ ■ (2) 



E(R,r],£) = E(R,ri,£ = 0)+£ a Eg\R,7i) 
^£ a £ /3 E%(R,rj) + ^ 
where we have introduced the compact notation 

bW ( r )= yg(R.i.£) 

Carried to all orders in £, this expansion is exact, but for practical purposes at sufficiently 
small fields we truncate this sum to define Ei(R, rj, £) as the sum of the first i + 1 terms in 
Eq. (g). We will also use similar expansions of R(£) and r\{£\ R and r] are the zero-field 
equilibrium structural parameters, routinely computed by using first-principles forces and 
stresses, and are assumed, along with the spontaneous polarization P s = P(Ro, f] , 0), to be 
known at the beginning of the tunability calculation. 

We first consider the collinear tunability of a polar insulator with a nonzero spontaneous 
polarization P s in a field £ parallel to P s , taken to define the z direction. This allows us to 
treat £, P, Xoo and x'oc as scalars, where Xoo(R, v) = E^(R, rj) and ^(R, rj) = E(%(R, rf). 
To clarify the logic of the presentation, we will write the expressions that follow for a single 
scalar R and a single component of the strain tensor rj; the generalization to the full set of 
atomic coordinates and strain tensor components is straightforward. 

We first assume that for any given structure, the DFPT package has the capability 
of calculating all third derivatives of the energy with respect to R, r\ and £. This leads 
us to consider E 3 and approximate x'oo(^v) ~ x'oo(-^o> Vo), resulting in the following key 
expressions: 

E 3 (R, V ,£) = E(R, V ,0) - P(R, V ,0)£ -l-Xoo(R,v)£ 2 -];xURo,Vo)£ 3 (3) 

2 o 

P{R^,£) = P(R, V ,0) + £ X oo(R,v) + lx'oo(Ro,Vo)£ 2 . (4) 



Our interest is in the £ derivative of the static dielectric susceptibility 



d Xs (£) d 2 P{£) d 2 P{R{£) 1 r 1 {£),£) 



d£ 



d£< 



d£< 



correct to leading order in the field. In the case under consideration, the zeroth order term 
suffices. This is a total derivative with respect to £, and thus involves derivatives of P with 
respect to R and rj as well as £ , as well as derivatives of R(£) and rj(£) at zero field, as 
follows: 
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For given £, the equilibrium structural parameters R{£ ) and rj(£) are found as solutions 
of the extremum equations: 








dE 3 (R, V ,£) dE(R, V ,0) dP(R,rj,0) £ 2 d X oo(R,v) 



dR 



dR 



-£ 



dR 



dR 



+ 



dE 3 (R,rj,£) dE(R, V ,0) g dPiR 1 r h 0) £ 2 d X oo(R,y) | 



dr\ 



drj 



drj 



drj 



(6) 
(7) 



We expand R{£) = R + R w £ + \R^£ 2 and r){£) = rj + + \4 2) ^ 2 - By solving the 
extremum equations to linear order in £, R^ and rf x * can be readily obtained in terms of 
second derivatives evaluated at the zero-field equilibrium structure: 
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At quadratic order, we similarly obtain R^ 2 > and rj( 2 \ 
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where 
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Inserting all these results into Eq. ^| completes the calculation of the linear dependence of 
Xs on electric field. 

Next, we assume that the highest derivatives of the energy that are available from the 
DFPT package used are the second derivatives ||. For example, we can directly compute 



|^|o but not |^r|o- As discussed in more detail below, a finite difference trick allows us 
to obtain relevant combinations of third derivatives one order higher in R or rj, but not in 



£ so that under these conditions, x'ooi^v)^ ^ ne electric field dependence of the electronic 
susceptibility, is not available. However, it is reasonable to neglect it, as the electronic 
susceptibility is small relative to the lattice contribution. To determine and r/ 2 ), we 
also need ^f-|o an d ^f^lo- From finite difference, we will be able to determine only the 
combination ^rlo-R^ 1 -* + ^Ho 7 / 1 ^- We obtain an estimate by assuming that all terms 
contribute equally to the sum. 

We first do a DFPT calculation at R , rjQ including all second derivatives. From this, 
we can obtain and 77 W from Eq. |[ Next, we chose a small field £\ and do a second 
DFPT calculation at structure R\ = R + R^£i and % = r q + r]^£ 1 . Using finite difference 
expressions, we can obtain the necessary combinations of third derivatives at the zero field 
equilibrium structure, as follows. For example, we need the £ dependence of Xoo(R(£), v{£)) 
only to linear order. The coefficient of the linear term is 

^| «W + ^l„-, (1 » (12) 

Since R\{£) — R$ = R^£ and — r]o= rf^S are linear in £ , we expand Xoo(R, v) around 
Ro, r/o to linear order in R — R and r\ — r] and evaluate at Ri, rjx. 

Xoo(R,v) = Xoo(Ro,vo) + ^-\oR (1 \R - Ro) + ^\oV {1) (v ~ Vo) 

to express the combination of derivatives given by Eq. |12] as 

(Xoo(Ri,rji) ~ Xoc(Ro,Vo))/£i) 

The calculation of P is a little more involved, since we need the £ dependence to quadratic 
order. 



8P, „m BP, 



P(R(£), V (£),0) = P(R , Vo , 0) + £(t^\ R {1) 
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The linear coefficient is a combination of available derivatives with R^ and r]^\ The two 
contributions to the quadratic coefficient (first derivatives of P and second derivatives of P) 
are determined by separate calculations. First, we use 

dP dP 
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to obtain 
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Then finally, we have to express and r/ 2 ) (from Eq. [|) in terms of available derivatives. 
This involves using finite difference expressions 
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Inserting all these results into Eq. [5] completes the calculation of the linear dependence of 
Xs on electric field. 



DISCUSSION 

The expressions presented in the previous section are for the collinear tunability of a ma- 
terial with a nonzero spontaneous polarization, where symmetry allows a linear dependence 
of Xs on £. To obtain tunability of paraelectrics and noncollinear tunability of ferroelectrics, 
where the linear term is zero by symmetry, requires going one order higher in the electric 
field. To this end, the necessary combinations of fourth-order energy derivatives can be 
obtained from DFPT calculations of third derivatives by an extension of the finite differ- 
ence approach described above to obtain the necessary combinations of third derivatives 
from calculated second derivatives. As in the latter calculation, the highest derivatives of 
the energy with respect to field cannot be obtained by finite difference calculations in the 
structure, and have to be estimated. Such approximations are acceptable in light of the fact 
that the lattice contribution to the dielectric constant dominates. 

Extension to higher order is also of interest for collinear tunability of ferroelectrics, as 
measurements in films allow the determination of the electric field dependence of Xs up to 
very high fields, showing strongly nonlinear behavior that cannot be investigated using 
the relatively low-order expansions of the present work. The development of methods for 
density functional calculations in finite fields |TD] will make it possible to investigate this 
regime. However, for the low-field tunability, the method described here should be practical 
for even fairly complex-structured materials. 
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